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FIBER-BASED OPTICAL LOW 
COHERENCE TOMOGRAPHY 

CROSS-REFERENCES TO RELATED APPLICATIONS 
[0001] This application claims the benefit of U.S. Provisional Patent Application 
No. 60/251,658, filed January 12, 2001. This application is a continuation-in-part of U.S. 
Patent Application 10/044,421, filed January 11, 2002, and incorporated herein by 
5 reference. 

BACKGROUND OF THE INVENTION 

[0002] The present invention relates generally to the field of optics and in particular to 
an interferometric system for birefringence characterization of one or more samples. 

[0003] Precise characterization of static or dynamic material birefringence is important 
1 0 and useful in applications that range from testing of optical components to measuring the 
viability and function of biologic tissue. While there are techniques that measure linear 
birefringence of samples, there remains a need for improved resolution and 
characterization, especially for samples or materials that are highly transparent, turbid 
and/or highly scattered. For example, fiber-based optical low-coherence reflectometery 
1 5 (OLCR) systems have been found to be rugged, relatively alignment-free, may be 
integrated with other sensing platforms and are flexible in interrogating a variety of 
samples. Unfortunately, fiber implementation of a polarization-sensitive OLCR system 
with standard single mode fiber is complicated since polarization state in a single mode 
fiber randomly fluctuates due to microscopic defects in the fiber, core ellipticity and 
20 external environmental perturbations. In addition, non-polarization maintaining (PM) 
fiber-based polarization-sensitive OLCR systems that perform birefringence 
characterization, work with multiple measurements and with different input polarization 
states to overcome the problem of fluctuating polarization state of light incident on the 
sample. Non-PM fiber-based setups cannot reliably measure birefringence of a sample 
25 with a single measurement even when using short lengths of fiber in the interferometer. 
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[0004] A technique that appears to be of more benefit for characterizing the 
birefringence of a sample is the polarization-sensitive optical low-coherence 
reflectometery (PS-OLCR). PS-OLCR is extremely useful for depth resolved 
birefringence characterization of transparent and highly scattering samples. With this 

5 technique, polarization maintaining single mode fibers can be used to construct a 
polarization-sensitive OLCR system, because the orthogonal polarization modes are 
isolated and preserve the input polarization state. Multiple measurements for 
birefringence characterization, however, limit detection of fast transient birefringence 
changes in a sample and are subject to any motion related artifact between multiple 

10 acquisitions. Although PM fiber based polarization-sensitive OCLR have been proposed 
to solve such problems, functioning and working systems have not been adequately 
constructed. Another disadvantage of current PM-OLCR systems is that they use bulk 
interferometric setups to interrogate samples for birefringence measurement. Therefore, 
there remains a need for a proper and efficient functioning PM-OLCR system useful for 

1 5 accurate birefringence characterization of a sample. 

SUMMARY OF THE INVENTION 

[0005] The present invention solves many problems associated with birefringence 
characterization, especially for samples with high scatter, turbid and/or transparent 
properties or of a sample in a solution with turbid or transparent properties. 

20 [0006] The present invention is based on the use of PM fiber-based dual channel OLCR 
to accurately measure birefringence, aspects of which have been disclosed by the 
inventors, disclosures of which are incorporated herein by reference. (See Dave DP and 
Milner TE. Electron Lett. 2001; 37:215; Dave DP and Milner TE, Opt Lett. 2000;25:227; 
both incorporated herein by reference.) In the previous disclosures, measured 

25 birefringence was sensitive to the rotation of the sample in the plane perpendicular to 
ranging. Here, configuration of a PM fiber-based dual channel OLCR system for 
birefringence characterization of one or more samples is described, wherein 
characterization may be made with a single measurement that is also insensitive to sample 
rotation. 

30 [0007] Generally, and in one form of the present invention, is a polarization- 
maintaining fiber-based polarization sensitive optical low coherence reflectometer for 
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depth resolved birefringence measurement. Unlike other fiber based polarization- 
sensitive optical low coherence reflectometers, linear birefringence of a sample can be 
measured from data recorded in a single A-Scan. Alternatively, multiple measurements 
may be recorded from a single sample and may provide substantially better measure of 
5 the specimen birefringence. The present invention provides for the simultaneous 
measurement of retardation and orientation of birefringent axes, where measured 
retardation is insensitive to sample rotation in the plane perpendicular to ranging. 

[0008] More specifically, one form of the present invention is a fiber-based optical low- 
coherence reflectometer comprising a polarization-maintaining source path, reference 

10 path, sample path, and detection path, wherein the polarization-maintaining sample path 
is optically aligned with a collimating lens, a variable wave retarder, and a focusing lens 
and the focusing lens is disposed to focus light on a sample. The polarization- 
maintaining source path, reference path, sample path and detection paths are each 
connected to a polarization-maintaining path coupler. The polarization-maintaining path 

1 5 coupler separates light into polarization-maintaining sample and reference paths while 
maintaining energy separation of optical signals. The polarization-maintaining detection 
path is optically aligned to provide for two mutually orthogonal light beams, each capable 
of producing an output signal that may be detected and analyzed to provide data (e.g., 
birefringence data) about the sample. The sample may include a turbid or transparent 

20 sample as well as a microfluidic chip. 

[0009] In another form, the invention provides for a method for characterizing 
birefringence of a sample comprising the steps of (1) creating a polarization-maintaining 
optical source path using a broadband light source; (2) creating a polarization-maintaining 
optical reference path that is optically coupled to a first collimator directed to an optical 
25 delay line with dispersion control; (3) creating a polarization-maintaining optical sample 
path that is optically coupled to a second collimator, a variable wave retarder, and a 
focusing lens, wherein the focusing lens focuses light on the sample; (4) creating a 
polarization-maintaining optical detection path optically coupled to a third collimator and 
a polarizing beam splitter, wherein the polarizing beam splitter is optically coupled to a 
30 first and second photodetectors that produce a first and second output signal, respectively; 
(5) connecting the polarization-maintaining source path, reference path, sample path and 
the detection path to a polarization-maintaining path coupler; (6) converting the first and 
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second output signals from the polarization-maintaining optical detection path with an 
analog-to-digital converter; and (7) connecting a processor to the analog-to-digital 
converter for collection of birefringent data about the sample. 

[0010] As with other forms of the invention, the apparatus and method of the present 
5 invention also provides for a polarization-maintaining optical fiber sample path optically 
aligned with a collimating lens, a variable wave or quarter wave retarder, and a focusing 
lens, wherein the focusing lens is disposed to focus light on a sample to characterize 
birefringence about the sample with rotation insensitivity of the measured retardation of 
the birefringent sample and wherein light back scattered from the birefringent sample 
10 after traversing through the wave retarder is elliptically polarized. The polarization- 
maintaining optical fiber sample path may further comprise an optical catheter probe used 
for imaging and/or be configured to interrogate a sensor. 

[00 1 1 ] In yet another form, the present invention provides for a system of 
characterizing birefringence of a sample comprising a broad bandwidth optical light 
1 5 source; a polarization-maintaining optical source path incorporating a polarizing element 
that correlates optical signals in fast and slow fiber polarization channels and optically 
connects both channels to a polarization-maintaining path coupler; a polarization- 
maintaining path coupler that separates light into polarization-maintaining optical sample 
and reference paths while maintaining energy separation of optical signals in the fast and 
20 slow fiber polarization channels; a polarization-maintaining optical reference path 

optically connected to the polarization-maintaining path coupler and optically coupled to 
an optical delay line; a polarization-maintaining optical sample path optically connected 
to the polarization-maintaining path coupler, wherein the polarization-maintaining optical 
sample path comprises a quarter wave retarder and a focusing lens, wherein the focusing 
25 lens is disposed to focus light on the sample, said sample placed in front of the 

polarization-maintaining optical sample path from which birefringence is characterized; a 
polarization-maintaining optical detection path optically connected to the polarization- 
maintaining path coupler and a polarizing beam splitter that is optically coupled to a first 
and second photodetectors that produce first and second output signals, respectively, 
30 wherein the first and second output signals are filtered and amplified, an analog-to-digital 
converter connected to the filter-amplifier and a processor connected to the analog-to- 



4 



Attorney Docket No. 1 19927-1068 



PATENT 



digital converter. Data collected from the processor about the sample is then used to 
characterizing sample birefringence. 

[0012] One of many advantages of the present invention is that it provides a mechanism 
by which to obtain additional contrast of a sample. For example, the present invention 

5 may be used to image and/or identify structural components of a sample or material, 
especially in materials (e.g., biologic materials) whereby functional information is 
associated with transient changes in birefringence. In addition, the present invention may 
be used to monitor, in real-time, surgical procedures, especially those that rely on 
photothermal injury mechanisms or changes in subsurface tissue components. The 

10 present invention may also be used as a feedback signal to control laser procedures in real 
time (e.g., dosimetry) and in catheters using optical fiber to deliver light. 

[00 1 3] Those skilled in the art will further appreciate the above-noted features and 
advantages of the invention together with other important aspects thereof upon reading 
the detailed description that follows in conjunction with the drawings. 

1 5 BRIEF DESCRIPTION OF SEVERAL VIEWS OF THE DRAWINGS 

[0014] For a more complete understanding of the features and advantages of the present 
invention, reference is now made to the detailed description of the invention along with 
the accompanying figures in which corresponding numerals in the different figures refer 
to corresponding parts and in which: 

20 FIGURE 1 depicts a schematic in accordance with one aspect of the present 

invention; 

FIGURE 2 depicts a schematic in accordance with another aspect of the present 
invention; 

FIGURE 3 depicts the demodulated interference intensity signals recorded from 
25 two surfaces of a sample at 1 300 nm in accordance with one aspect of the present 
invention; 

FIGURE 4 depicts the measured single-pass phase retardation versus the angular 
rotation of a sample in accordance with one aspect of the present invention; 
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FIGURE 5 depicts the measured and actual orientation of the fast axis of a sample 
with respect to its horizontal axis in accordance with one aspect of the present invention; 

FIGURE 6 depicts an example of a non-linear algorithm to determine depth- 
resolved phase retardation in accordance with one aspect of the present invention; and 

5 FIGURE 7 depicts one aspect of a procedure to determine # z , A z) using as input, 

<p f ,e f ,e 0 estimates determined from the Levenberg-Marquardt non-linear algorithm and 
recorded PS-OCT data. 

DETAILED DESCRIPTION 

[00 1 5] Although making and using various embodiments of the present invention are 
10 discussed in detail below, it should be appreciated that the present invention provides 
many inventive concepts that may be embodied in a wide variety of contexts. The 
specific aspects and embodiments discussed herein are merely illustrative of ways to 
make and use the invention, and do not limit the scope of the invention. 

[001 6] In the description which follows like parts may be marked throughout the 
1 5 specification and drawing with the same reference numerals, respectively. The drawing 
figures are not necessarily to scale and certain features may be shown exaggerated in 
scale or in somewhat generalized or schematic form in the interest of clarity and 
conciseness. 

[0017] FIGURE 1 depicts one aspect of the present invention, the PM fiber 
20 polarization-sensitive OLCR system, also referred to as the PM fiber-based polarization 
sensitive optical coherence tomography (PS-OCT). A light source 10 emits broadband 
light that is polarized by a polarizer 12 and coupled to the fast or slow axis of an input 
port fiber of a path coupler 14. In one embodiment, the light source includes an optical 
semiconductor amplifier centered at 3^=1 .3 1 um with a coherence length of 24 urn. Here, 
25 the polarizer 12 is a fiber bench polarizer that polarizes the emitted light followed by its 
being coupled into the fast or slow axis of the input port fiber of a 2 x 2 PM coupler. A 
polarizing extension fiber spliced at 0° to the input fiber of the 2 x 2 PM coupler may be 
utilized to provide a high extinction ratio between light in orthogonal modes of the PM 
fiber 
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[001 8] In general, the fast and slow axes of the input port fiber form the two OLCR 
polarization channels, Channel 1 (Chi) and Channel 2 (Ch2), respectively. Back 
reflected light from a reference path and sample path in each channel mix at the 2 x 2 
coupler to form interference signals. 

[0019] In a PM fiber interferometer, interference signals formed in Chi and Ch2 will 
overlap only if optical path-length differences between sample and reference paths for 
each channel are equal. As shown in FIGURE 1, the length of the coupler fiber 16 in the 
reference path and length of the coupler fiber 18 in the sample path are equal. In 
addition, light in Chi and Ch2 travel the same optical path-length in the reference path 
after a 45° splice. As such, fibers attached to the coupler fiber 1 8 in the sample path 
ensure that light in Chi and Ch2 propagate equivalent optical path-lengths. 

[0020] Owing to the birefringence of the PM fiber, any single segment of PM fiber 
attached to the coupler fiber in the sample path will result in unequal optical path-lengths 
for the two channels. Now referring to the sample path in FIGURE 1 , when two equal 
length segments of PM fiber are spliced at 90° to each other, total optical path-length 
traversed by light in Chi and Ch2 is equal. In one embodiment of the sample path, there 
are two segments spliced to the PM fiber: segment A and segment B. One end of 
segment A 20 is spliced at 0° to the coupler fiber 1 8 and the other end is spliced to 
segment B 22 at 90°. In general, the length of segment A 20 is equivalent to that of 
segment B 22. Here, a 90° splice interchanges the two channels, i.e., Chi completely 
couples into the Ch2 and vice-versa and provide equal path-length between channels. 
Segment B 22 is terminated with a connector 22 that may be an angle-cleaved connector. 
The connector 22 is attached to a collimator 24. 

[0021] The sample path also includes a quarter wave or variable wave retarder 38, such 
as a achromatic quarter wave retarder, a liquid crystal variable retarder (LCVR), or 
possibly a short segment of PM fiber spliced at 45 degrees to the main fiber with a 
piezoelectric stress element to vary the birefringence. In one embodiment, the quarter 
wave retarder is an achromatic quarter wave plate. Importantly, the quarter wave retarder 
provides a quarter wave of retardation to each of the OLCR source wavelengths. The 
quarter wave retarder or variable retarder 38 is generally aligned with its slow axis 
oriented at 45° to the fast or slow axis of the PM fiber so that circularly polarized light is 
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incident on the sample. With use of a variable wave retarder, wave variation is generally 
is from zero to one wavelength. 

[0022] A focusing lens 58 focuses the polarized light to a single point on a sample. 
Light back scattered from the birefringent sample after traversing through the variable 
5 wave retarder 3 8 is elliptically polarized. 

[0023] Now referring to the reference path in FIGURE 1, output of the coupler fiber 16 
is spliced at 45° to the input fiber 26 of a phase modulator 28. In one embodiment, the 
input fiber is a PM fiber and the phase modulator 28 is a Lithium Niobate (LiNb0 3 ) 
waveguide electro-optic phase modulator. With this modulator, only light propagating 
10 along Ch2 of the input fiber 26 pigtailed to the phase modulator 28 is launched into the 
TE mode of the integrated optic waveguide. An output fiber 30 (e.g., PM fiber) from 
phase modulator 28 is spliced at 0° to a fiber segment 32 (e.g., PM fiber) that terminates 
with a connector 34 attached to a fiber collimator 36. The collimator 36 collimates light 
emitting from the optical fiber into an optical delay line with dispersion compensation. 

1 5 [0024] With the present invention, collimated light from the output fiber 30 may enter 
an optical delay line that provides path-length modulation and compensates for dispersion 
imbalance in the interferometer caused by the phase modulator 28. For example, the 
dispersion compensating delay line may consist of a pair of diffraction gratings 54 and a 
movable end mirror 56. The function of the grating pair is to compensate chromatic and 

20 waveguide dispersion introduced by the phase modulator 28. If left uncompensated, 
dispersion increases the width of the coherence function tenfold thereby drastically 
reducing the longitudinal resolution of the interferometer. After a round trip, input light 
traveling along one fiber axis in the reference path is equally coupled to both fiber axes at 
the 45° splice. 

25 [0025] Now referring to the detection path of FIGURE 1, interfering light that exits the 
fiber is collimated by collimator 43 and split into two linearly polarized light beams or 
paths 40 and 42 that are mutually orthogonal from each channel by a polarizing beam 
splitter 44. In one embodiment , the polarizing beamsplitter is a Wollaston prism that 
separates the two signal channels. Each linearly polarized light path is focused into and 

30 detected by a photodetector; polarized light path 40 is detected by photo-detector 46 and 
polarized light path 42 is detected by photodetector 48. Photo-detectors 46 and 48 are 
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each connected to an analog-to-digital converter (e.g., two-channel 12 bit) and signal 
output from each detector may be stored in computer memory for processing. The phase 
of each polarization channel may be determined by computing the arctangent of the 
filtered signal divided by its Hilbert transform. 

5 [0026] Still referring to FIGURE 1 , in the detection path, interfering light is split into 
two mutually orthogonal linearly polarized beams from each channel by a Wollaston 
prism 44 (e.g., large-angle Wollaston prism) and is detected by two photodetectors 46 and 
48. If the sample is a linearly birefringent material with retardation 8(z) and its fast axis 
is oriented at an angle 6 from the horizontal, the interference intensity signals detected in 
10 Chi and Ch2 can be written shown below: 

l l (z) = k iy [Rjz)cxp[-(Az/l c ) 2 ] 

xsin(5)cos(4*Az//l o + 20) (4) 

15 I 2 (z) = * 2 V^)exp[-(Az/0 2 

x cos(<5)cos(4^Az/A 0 ( 5 ) 
where Az is the optical path-length mismatch between the reference and the sample paths 

20 l c is the coherence length of the source, R s is the sample' s reflectivity, and k, and k 2 are 
proportionality constants that depend on the quantum efficiency of the detector and on the 
efficiency of light coupling into the detectors. In deriving equations (4) and (5), the light 
source's power spectrum is generally assumed to be Gaussian. In addition, it is assumed 
that the relative shift of the coherence function between two polarizations that is due to 

25 the sample's birefringence is negligible, and that the orientation of the axes remains fixed 
as a function of sample depth. 

[0027] The phase function of the interference intensity signal in equations (4) and (5) is 
given by the following equation: 

* m (z) = arctan[//[/ m (z)]/ m (z)}, m = 1,2 (6) 

30 

and the amplitude function is given by: 

A m (z) = ({H[I m (z)]r + [l a (z)} 2 ) in W 

3 5 [0028] An important feature of the present invention is the ability to make a 

measurement of specimen birefringence using only one incident state. In addition, 
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measurements using multiple incident polarization states may be made and will provide a 
substantially better measure of the specimen birefringence. To measure the birefringence 
using multiple incident polarization states, the present invention includes the following: 
(a) when using a quarter wave retarder in the sample path, it should be replaced with a 
variable wave retarder (e.g., a liquid crystal variable retarder or a short segment of PM 
fiber with external stress elements); and (b) analysis of birefringence may be generalized 
to the Poincare sphere. Examples of this are provided in one of the following sections. 

[0029] Now referring to FIGURE 2, the present embodiment provides a fiber-based 
optical low-coherence reflectometer 100 comprising a path coupler 1 10 that separates 
light into a birefringent sample 140 and a reference path 130 while maintaining energy 
separation of optical signals in the fast and slow fiber polarization channels. As used 
herein, connections, alignments, optical connections and optical couplings may be used 
interchangeably. The reflectometer includes a source path 120 comprising a first 
polarization-maintaining optical fiber 122 having a first end and a second end, wherein 
the first end of the first polarization-maintaining optical fiber is optically coupled to a 
light source 124 and the second end is connected to a polarizer 126 that splits the light 
source into a first and second polarization channels with independent phase components. 
The polarizer 126 is connected to a second polarization-maintaining optical fiber 128 
having a first end and a second end, the first end connected to the polarizer 126 and the 
second end connected to the path coupler 1 10. In one aspect, the polarizer output fiber is 
pigtailed to a polarization-maintaining optical fiber spliced at 0° with termination at the 
path coupler 110. 

[0030] Still referring to FIGURE 2, the reference path 1 30 comprises a third 
polarization-maintaining optical fiber 132 and fourth polarization-maintaining optical 
fiber 134. The third polarization-maintaining optical fiber 132 has a first end and a 
second end, the first end connected to the path coupler 1 10, the second end connected to a 
phase modulator 136. The fourth polarization-maintaining optical fiber 140 having a first 
end and a second end, the first end connected to the phase modulator 136 and the second 
end to a connector 137 and optically aligned with a first collimator 138 that collimates 
light emitting from the second end of the fourth optical fiber into an optical delay line 
139. In one embodiment, the output fiber from path coupler 1 10 is spliced at 45° to the 
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input fiber of phase modulator 136. In addition, the fourth polarization-maintaining 
optical fiber 134 is spliced at 0° before entry to connector 137. 

[003 1] Still referring to FIGURE 2, the sample path 140 comprises a fifth polarization- 
maintaining optical fiber fiber 142 having a first and a second end, the first end connected 
to path coupler 1 10, the second end to a connector 144 and optically aligned with a 
second collimator 145 that collimates light emitting from the fifth optical fiber 142 to a 
variable or quarter wave retarder 146 that ends on sample 149. Sample path 140 further 
comprises two segments, segment A 147 and segment B 148. One end of segment A 147 
is spliced at 0° to the path coupler 1 10 output fiber and the other end is spliced to segment 
B 148 at 90°. Segment B is terminated with connector 144. Fiber lengths of segments A 
147 and B 148 may be equivalent. The 90° splice interchanges the two channels. 

[0032] Still referring to FIGURE 2, the invention also includes a detection path 
comprising a sixth polarization-maintaining optical fiber 152 having a first end and a 
second end, the first end connected to the path coupler 110, the second end aligned with a 
third collimator 154 that collimates light emitting from the sixth optical fiber 152 onto a 
polarizing beam splitter 156. The polarizing beam splitter 1 56, such as a Wollaston 
prism, splits the light from the sixth polarization-maintaining optical fiber 152 into a first 
beam 158 and a second beam 160 that are mutually orthogonal and capable of producing 
a first and second output signal. 

[0033] In another aspect of the present invention, the first output signal is provided by a 
first photodetector that detects light from the first beam; the second output signal is 
provided by a second photodetector that detects light from the second beam. As such, a 
first output signal pass from the first photodetector and a second output signal pass from 
the second photodetector, each pass through a bandpass filter and amplifier to produce a 
first and a second filtered signal, respectively. The two photodetectors are connected to 
an analog-to-digital converter (e.g., two-channel 12 bit) and the first and second filtered 
signals from the analog-to-digital converter are analyzed immediately or stored in a 
processor (e.g., as computer memory) for later processing. Frequency multiplexed 
signals may be separated by digital filtering and the phase of each frequency channel is 
calculated by computing the arctangent of the filtered signal divided by its Hilbert 
transform. 
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Examples Using the PS-OLCR System 

[0034] In one example, the PS-OLCR system uses a mica quarter wave plate (QWP) set 
at 1300 nm and placed as a sample. Collimated light in the sample path is focused with a 
40 mm focal length achromatic doublet lens on to the sample mica QWP. Scanning in 
depth, interference intensity signal from the front and back surface of the 100 um mica 
QWP in Chi and Ch2 are recorded. Demodulated amplitudes of interference intensity 
signals from the front and back surface of the mica QWP in each channel are shown in 
FIGURE 3. As expected, going from front to back surface of mica QWP, the signal 
completely shifts from Chi to Ch2. Due to the non-ideal nature of the QWPs and finite 
cross-coupling isolation between polarization channels in the PM fiber, residual signal is 
present in Chi and Ch2 from light reflecting from front and back surface of the mica 
QWP. 

[0035] To demonstrate rotation insensitivity of the measured retardation of a 
birefringent sample a QWP at 850 nm (8 « 58° at 1300 nm) mounted on a rotation stage 
was placed in the sample path. Single-pass retardation from passage through the sample 
for various angular orientations was calculated. FIGURE 4 shows measured phase 
retardation versus sample rotation from O 0 to 180°, plotted in increments of 5°. Here, the 
average value of the measured phase retardation over the whole rotation range was 58.4° 
with a standard deviation of 1 .2°. 

[0036] Calculating the phase difference between the two channels gives the orientation 
(20) of the slow axis of the sample with respect to the horizontal. Shown in FIGURE 5 
are the measured and expected values of 2d as a function of QWP rotation. Since the PM 
fiber segments in the sample and reference paths are not perfectly matched, an offset 
phase difference exists between Chi and Ch2. The measured values of 9 plotted in 
FIGURE 5 have been corrected for the offset phase difference. 

Non-linear Algorithm to Determine Depth-Resolved Phase Retardation 

[0037] With the present invention, optical elements in the reference and sample paths 
systematically alter horizontal and vertical interference fringe intensities. By use of the 
present invention, these altered effects may be readily calibrated and corrected to 
determine the true polarization state of light backscattered from a sample. For example, 
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with the reference path, mismatch of fiber length (between reference and sample paths) 
modify the relative amplitude (r=A^A H ) and offset the relative phase (fr=fa-v) of 
reflected light. Likewise, in the sample path, a variable wave retarder, such as a liquid 
crystal variable retarder, offsets the relative phase (fa-v) of transmitted light by the 
5 applied retardation S LC vr • Subsequence relative phase effects are corrected by 
multiplying the analytic signal for the horizontal interference fringe intensity by a 

correction factor re'^ _<W . 

[0038] Further determination of the depth-resolved phase retardation [A%z,Az)] in a 
form-birefringent sample requires a non-linear algorithm that uses recorded data about a 

1 0 sample as input to estimate the pseudo fast axis orientation F(<fy, 0/) and cone apex-angle 
(0 O ). Determination of 0 f ) and 6 0 is problematical because recorded data about the 
sample is degraded by speckle noise. Moreover, because recorded data about the sample 
is naturally interpreted on the Poincare sphere, underlying geometry of the problem to 
determine F(fo 6j) and 6 0 is spherical and a nonlinear algorithm must be applied to 

15 determine AS(z,Az). 

[0039] As such, the present invention provides a non-linear algorithm that consists of 
three components (FIGURE 6): (A) initial estimate of pseudo fast axis orientation 
[$O=0),0/i=0)] and cone apex-angle [9&=0)\, (B) determination of F and 9 0 \ and (C) 
computation of the depth-resolved phase retardation [S(z,Az)]. 

20 Examples of Pseudo Fast Axis Orientation and Cone Apex- Angle Estimates 

[0040] Because non-linear algorithms can give erroneous results when local minima 
exist in the residual function, determination of good initial estimates for the pseudo fast 
axis orientation [$(/=0), 0/0=0)] and cone apex-angle [ft(i=0)] are important for 
convergence. An approach to determine [<p f (i = 0),0 f {i = 0),6 o (.i = 0)] directly from 
25 recorded data about a sample must manage the degrading effect of speckle noise. 

[004 1 ] The present invention provides a simple approach to determine an initial 
estimate for [F(^0/),, 0J while managing and suppressing effect of speckle noise. 
Normalized Stake's parameters (S Q , Su, S v ) determined from recorded data about a 
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sample are plotted versus depth z and fitted to a polynomial. Degree (2 nd or 3 rd ) of the 
fitting polynomial is selected based on prior knowledge of the phase retardation (&«)• 

[0042] Polynomials that best fit (S Q , S Ut S v ) give the fitted Stoke's parameters and are 
denoted by [Soj(z\SijJ(z),SvJ(z)]. As used herein, phase retardation of the birefringent 

5 layer is 5 es , and uses the fitted Stokes parameters [Sq^z^Su^JSv/z)] to compute 

[<p f (i = O),0 f (i = 0),& o (i = 0) ]. First, the fitted Stoke's parameters [S&(z),Suj(z)£vj(z)] are 
divided into N s equal-length (Az/N s ) intervals where Az is optical thickness of the form- 
birefringent sample. Second, a representative Stoke's vector (Si) for each interval is 
computed from the fitted parameters [Sq{z),Su{z),Sv{z)]. Third, using 8m and a pair of 

1 0 representative Stoke's vectors (Si and S m ) corresponding to intervals / and m, a fast axis 
(n, m ) and cone apex-angle (9 0 ,i m ) are uniquely determined. By considering all pairs of 
representative Stokes vectors (Si and S m ), initial estimates of the fast axis (<«>) and cone 
apex-angle [0^i=O)] are formed by computing the weighted mean over n tm and 0 o ,i m 
respectively (Equation 8). The weighting factor (w /m ) is given by a simple power of the 

1 5 relative distance (e.g., \l-m\ i/2 ) between / and m 'th intervals. 

< « >= { I w lm n lm 1/ I w lm , 6 0 (i = 0) = | I w lm -e 0tlm \l I w lm (g) 

\l*m J l±m \l*m J l*m 

[0043] Initial estimates <p f (i = 0) and 6 f (i = 0) corresponding to <n> are determined 
directly using the representation ( («) ^sin^cos^Sg+sin^sin^+cos^) Q f a unit 
20 vector in spherical coordinates. 

Determination of Pseudo Fast-Axis Orientation and Cone Apex- Angle 

[0044] For a statistically unbiased estimate, the pseudo fast axis orientation F(<fe Of) and 
cone apex-angle (0 O ) must be determined using recorded data about a sample without the 
bias introduced by fitting. 

25 [0045] The present invention provides a non-linear algorithm that determines an 
unbiased estimate of [F(<fr, 6j), 0 O ] by minimizing a residual function (R; Equation 9) 
formed by the composite sum of distances (e,) on the Poincare sphere formed between the 
PS-OCT data points (Si) and the arc specified by [<fr, Of, 0„]. Here, s is the shortest 



14 



Attorney Docket No. 1 19927-1068 



PATENT 



10 



15 



20 



distance between the i'th PS-OCT data point (Si) and the arc specified by [fa, 0 f , 0 O ]. The 
distance e t is given by the angle along the great circle that is perpendicular to the arc and 
contains S,\ 



[0046] The pseudo fast-axis n(fa, 0j) and cone apex-angle (0 O ) are determined by 
minimizing the residual function [R(<p f ,0 f ,0 o )] using a Levenberg-Marquardt nonlinear 
algorithm. The algorithm starts with the initial estimate <p f (i = O),0 f (i = O),0 o (i = 0) derived 
from the fitted Stokes parameters (see previous section) and uses an iterative approach to 
determine q>f,0f,0 o . 

Computation of Depth-Resolved Phase Retardation 

[0047] The present invention provides a procedure to determine d(z,Az) using as input, 
<p f ,0 f ,0 o estimates determined from the Levenberg-Marquardt non-linear algorithm and 
recorded data about a sample (Si). In this procedure, a plane is formed that is both tangent 
to the Poincare sphere and perpendicular to the pseudo fast axis n(fa, 0$. Recorded PS- 
OCT data points (Si) are projected onto the tangent plane along the line extending from 
the center of the sphere to S\ (FIGURE 7 A). A right-handed orthogonal coordinate 
system (u,v) is defined on the tangent plane with an origin corresponding to the pseudo 
fast axis [F (fa 0/)] (FIGURE 7B). The angle (ai) formed between each projected data 
point (S\) and the u axis in the tangent plane is computed and plotted against the data 
point index L Slope (m) of the linear least squares fit line to the (i,ai) data pairs gives the 
phase retardation per unit index change (FIGURE 7C). 

[0048] With the present invention, the depth-resolved phase retardation [S(z,Az)] is 
given by d(z,Az)=N p m where N p is the number of data points recorded over optical depth 
Az. 

[0049] Thus, the present invention is capable of depth resolved birefringence 
characterization with a single and with multiple measurements. The design of the 
reflectometer and optical coherence tomography system of the present invention places 
no limit on the length of fiber in the sample path as long as it is matched with proper 
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length of PM fiber in the reference path. Inclusion of the electro-optic modulator not only 
provides a stable carrier frequency but also permits measurement of fast transient 
birefringence. Cross-coupling between PM fiber axes, either due to finite isolation or 
imperfect splice angle between fiber axes, may limit the dynamic range and birefringence 
5 measurement accuracy of the PM fiber based PS-OLCR instrument. 

[0050] There are several advantages of the present invention that include: (a) an 
alignment-free system; (b) ease of use; (c) ability to integrate with other systems (e.g., 
microchips, sensors, catheters, lasers); and (c) an ability to probe sites that cannot be 
interrogated with current bulk PS-OLCR systems, such as for birefringence imaging of 

10 locations in a mammalian body. For example, the present invention may be used to probe 
a microfluidic chip placed in the sample path. The sample path may also consist of an 
optical catheter probe for imaging a sample, such a mechanical, electrical or biologic 
materials (e.g., an organ, tissue, or luminal wall of a coronary artery). The sample path 
may be configured to "read" or "interrogate" a sensor, such as a surface plasmon 

15 resonance device. 

[005 1] Additional objects, advantages and novel features of the invention as set forth in 
the description, will be apparent to one skilled in the art after reading the foregoing 
detailed description or may be learned by practice of the invention. The objects and 
advantages of the invention may be realized and attained by means of the instruments and 
20 combinations particularly pointed out here. 



